Mobile DNA

Bacterial insertion sequences (IS) are compact transposable elements that encode proteins required for their mobility and maintenance, yet many also encode accessory proteins with poorly understood functions. For example, IS605-family elements often encode a transposase called TnpA and an RNA-guided nuclease called TnpB that supports transposon maintenance, alongside an additional ribbon-helix-helix protein named Arc. Though the roles of TnpA and TnpB have been extensively studied in recent years, the enigmatic function of Arc has not been investigated. Here, we show that Arc acts as a transcriptional repressor to directly bind the transposon's native promoter sequence regulating TnpA and TnpB gene expression. By systematically testing Arc-containing IS605 elements, we identified a conserved binding pattern at intergenic transposon sequences neighboring protein-coding genes through chromatin immunoprecipitation and sequencing analyses. We then used fluorescence reporter assays and demonstrated that these intergenic sequences function as strong promoters, and that the presence of Arc dramatically reduces their gene expression. Together, these findings identify Arc as a transposon-encoded transcriptional repressor, revealing a regulatory layer that may promote long-term persistence of IS605-family elements by keeping their activity in check. The widespread association of Arc homologs with diverse mobile elements and cellular genes suggests that these compact regulators may more broadly restrain the expression of neighboring genetic machinery across varied genomic contexts. neighboring genetic machinery across varied genomic contexts.

Transposable elements (TEs) are a major source of genomic diversity in mammals, yet their regulatory roles in the bovine genome remain poorly understood. Through characterizing bovine TE landscape, despite the substantial proportion (25.6%) of ruminant-specific TEs, we observe age- and class-dependent genomic distribution patterns similar to those observed in other mammals. Next, we profile TE and gene expression dynamics in pre-implantation embryos generated in vivo (IVV), by in vitro fertilization (IVT), and through somatic cell nuclear transfer (SCNT). The zygotic genome activation (ZGA) is shifted from the 4-cell stage to the 8-cell stage in IVT and SCNT embryos compared to IVV embryos. SCNT embryos exhibit impaired initiation of early transcription programs at the 4-cell stage and disrupted developmental trajectories, including abnormal activation of pluripotency-associated genes. A subset of retroviral LTR elements are strongly activated at ZGA in IVV and IVT embryos, whereas their activation is markedly muted in SCNT embryos, suggesting that impaired gene and TE reprogramming may contribute to the developmental defects commonly observed in SCNT embryos. By epigenomic profiling, the MLT1D elements from the ERVL-MaLR LTR family lose repressive marks and gain H3K27ac at ZGA, together with DUXA-binding motif enrichment. Knockdown of DUXA in bovine embryos significantly reduced MLT1D expression and ZGA marker genes. We propose that a subset of DUXA-enriched MLT1D functions as enhancers that promote ZGA. Overall, our study provides new insights into the regulatory roles of TEs during bovine embryogenesis and establishes a framework for comparative studies of TE-mediated gene regulation in early mammalian development.

We are interested in DNA stiffness in vitro and in vivo. Classical approaches for measuring this parameter in vitro involve tedious T4 DNA ligase-mediated cyclization kinetics experiments performed for one DNA fragment length at a time to establish effective local end-end concentrations (J-factors). Protein-mediated DNA looping is a more relevant biological reaction. An efficient approach to determine length-dependent protein-mediated DNA J-loop values in vitro has been lacking. We describe LOOP-TAG, a massively-parallel method. We leverage intramolecular DNA looping and tagmentation catalyzed by Tn5 transposase (Tnp) tethered to the terminus of a ~1,000-bp bead-bound duplex DNA. We detected loops terminating with Tnp. Counts of the resulting DNA looping-dependent bead-bound tagmentation products from deep sequencing determine DNA loop length probabilities at single-base pair resolution for all bead-bound fragments. When normalized to DNA fragment length probabilities from tagmentation using known concentrations of free Tnp, length-dependent J-loop values are directly determined. We demonstrate LOOP-TAG for Tnp tethered in two ways and in the presence of DNA-kinking architectural DNA binding protein Nhp6A. Intermediate length J-loop values are consistent with expectations of wormlike chain theory after accounting for strong sequence-dependence of Tnp. DNA loops are shifted to smaller sizes in the presence of Nhp6A, as predicted.

Recombinant vectors derived from adeno-associated viruses (rAAVs) are a mainstay of human gene therapy. rAAVs are produced from plasmids containing transgene cassettes flanked by inverted terminal repeats (ITRs), which form structured DNA elements that stabilize the ends of the single-stranded viral genome and are the only viral sequences required in cis for genome packaging. For decades, it has been recognized that propagation of ITR-containing plasmids can result in deletions and other mutations, prompting the use of specialized bacterial strains, modified growth conditions, and truncated or altered ITRs. Despite these practices, ITR instability remains a persistent source of plasmid heterogeneity. To identify determinants of ITR stability, we evaluated ITR integrity in one of the original cloned AAV2 genome isolates, a reconstructed AAV2 plasmid, and a synthetic rAAV vector containing full-length native AAV2 ITRs. We established a quantitative bioinformatic workflow for analyzing ITR-containing plasmids and virus preparations from raw Oxford Nanopore sequencing data. These experiments showed that ITRs were highly stable during short-term culture, whereas prolonged culture revealed strong positional effects, with preferential loss or mutation of the ITR nearest the plasmid origin of replication. Consistent with this model, a survey of 7,041 sequence-verifiable AAV plasmids from the Addgene repository identified a widely disseminated 11-bp ITR deletion in 4,773 plasmids; among analyzable two-ITR plasmids, this deletion was located in the origin-proximal ITR in 95.3% of cases. Guided by these findings, we constructed a novel rAAV entry vector with stable full-length native AAV2 ITRs that enabled efficient packaging of a 4,750-bp all-in-one CRISPR-Cas9 cassette. Finally, we developed a cell-based strategy to compare the effects of ITR mutations on rAAV genome integration, providing preliminary evidence that ITR sequence variation can influence integration outcomes. Together, these findings show that ITR instability is a preventable, position-dependent property of plasmid architecture and identify ITR integrity as an important variable in rAAV vector design and quality control.

Insertion sequences (IS) are the most common type of transposable element in prokaryotes and shape the structure of genomes through transposition and by providing a substrate for recombination. Despite the mutational impact of IS, the evolutionary dynamics of most elements in host species remain unknown. Here we study the dynamics of IS6110 in 10,000 strains of the Mycobacterium tuberculosis complex (MTBC). We developed a tool that allows the detection and comparison of IS insertions from short reads without using a reference genome. Using ancestral state reconstruction (ASR) on presence-absence patterns of IS6110, we describe the distribution of copy numbers (CNs) in the MTBC, infer birth rates of the element, and identify genomic regions with large numbers of parallel IS6110 insertions. Copy numbers in the MTBC range from 1 in some clades to more than 30 in strains of La3 (M. orygis). IS6110 birth rates scale approximately linearly with copy number and are elevated on terminal branches, consistent with the delayed action of purifying selection. A key characteristic of IS6110 is its occurrence in hotspots: the 5% most frequently targeted regions account for half of all independent insertion events. The motif 5'-TCTCAAAW-3' is enriched around target sites and in hotspots, suggesting that the accumulation of insertions in these regions results through a combination of non-random insertion and purifying selection in other regions. To conclude the study, we propose a niche constraints model according to which the distribution of IS6110 in the MTBC is governed by the rarity of regions that have both suitable DNA properties and little functional value for the host.

Reliable classification of transposable elements (TEs) is important for genome annotation and evolutionary analysis, yet remains challenging: TE families differ widely across taxa, reference databases label them inconsistently, and many TEs are represented only by degraded copies. We introduce GLARE (Genomic Language model for Automated Repetitive Element classification), a TE classifier built by fine-tuning a genomic language model on harmonized PanTEon and Repbase reference data to predict TE superfamily across 32 superfamilies in 11 orders. Among the genomic language models compared, Nucleotide Transformer v3 provided the strongest representation. Built on this backbone, GLARE exceeded all nine classifiers in the PanTEon benchmark, reaching 0.966 weighted superfamily F1 on its 30-superfamily task. GLARE also showed stronger generalization to clades absent from training, consistently exceeding state-of-the-art classifiers on holdouts spanning three taxonomic ranks: the order Poales, the phylum Arthropoda, and the kingdom Fungi. Training on both PanTEon and Repbase reduced the performance drop seen when a single-source model was applied to the other database. Under simulated sequence decay, GLARE was the most robust classifier tested, and it performed best on independently curated degraded libraries. These findings show that genomic language models fine-tuned on harmonized multi-source references can classify TEs robustly under the taxonomic, database, and sequence-quality shifts that matter in practice. GLARE is available at https://code.jgi.doe.gov/NHaas/glare.

Long non-coding RNAs (lncRNAs) are a class of regulatory RNAs with lengths exceeding 200 nucleotides. The lncRNA PURPL (p53 upregulated regulator of p53 levels) is a p53-regulated gene. PURPL is an oncogenic lncRNA that represses p53 function; it is upregulated in many tumors and promotes cancer cell proliferation. Kaposi sarcoma-associated herpesvirus (KSHV) is the causative agent of all forms of Kaposi sarcoma (KS) and is tightly associated with primary effusion lymphoma (PEL). We found upregulation of PURPL in both de novo and chronically KSHV-infected cells. We found that p73 (a p53 homolog) associates with the PURPL promoter and upregulates PURPL expression. Previously, it has been shown that KSHV-encoded latency-associated nuclear antigen (LANA) stabilizes p73. We found that ectopic expression of LANA increased both p73 protein levels and PURPL expression. p73 knockdown abrogates PURPL upregulation following KSHV infection. Antisense oligonucleotides (ASO) targeting PURPL reduced PEL cell proliferation and increased acetylated p53 levels. In contrast to p53, p73 is not repressed by PURPL RNA. In summary, our study reveals how an oncogenic virus can utilize p73 to upregulate the oncogenic lncRNA PURPL as an additional way to repress p53. Moreover, we propose another mechanism by which p73 may promote tumor growth.

LINE-1 (L1) is transcribed by RNA polymerase (Pol) II, whereas the non-coding Short INterspersed Element (SINE) Alu is transcribed by Pol III and exploits L1 ORF2p for mobilization. Conventional L1 retrotransposition reporters rely on spliceosome-dependent introns and cannot directly monitor Alu retrotransposition. Here, we developed EGFPTet, a fluorescent SINE reporter in which EGFP is interrupted by a Tetrahymena thermophila group I self-splicing intron. EGFPTet enabled flow cytometry, fluorescence microscopy, cell sorting, and time-course analysis of Alu mobilization without colony formation. Reporter activation required ORF2p reverse transcriptase activity, and recovered Alu-EGFPTet insertions displayed retrotransposition hallmarks. EGFPTet quantified Alu retrotransposition across cell lines and host-factor perturbations and was repurposed to detect mobilization of the phylogenetically distinct mouse B2 SINE. Defined amino acid substitutions generated an mTurquoise2 derivative without redesigning the self-splicing junctions, allowing dual-color analysis of Alu trans-mobilization and L1 cis retrotransposition in individual cells. In this assay, wild-type L1 supported a higher Alu retrotransposition frequency than the ORF1p-mutant L1. Despite their shared dependence on ORF2p, Alu co-expression did not significantly reduce L1 retrotransposition, and Alu/L1 dual-positive counts exceeded independence-based expectations, indicating preferential co-occurrence of the two reporter signals. EGFPTet provides a fluorescent platform for investigating the shared and distinct cellular conditions that support SINE and L1 mobilization.

The 5'-untranslated region (5'-UTR) of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is highly structured and comprises five conserved stem-loops that serve as cis-acting regulatory elements. During infection, host proteins are recruited to these elements and regulate multiple steps during the viral life cycle. Stem-loop 5 (SL5) is highly conserved among different alpha- and betacoronaviruses and includes the start codon of non-structural protein 1 (nsp1). However, the function of SL5 in viral replication and its interacting proteins remain to be determined. Here, we identified heterogeneous nuclear ribonucleoprotein M (HNRNPM) as a host protein functionally interacting with SL5. HNRNPM translocates into the cytoplasm during SARS-CoV-2 infection and promotes protein synthesis from the viral 5'-UTR. Mutations in the lower stem of SL5 impair HNRNPM binding and reduce protein synthesis. Notably, the interaction between HNRNPM and SL5 is conserved among different alpha- and betacoronaviruses. Our findings indicate that HNRNPM is an important host factor that regulates viral translation to promote SARS-CoV-2 replication.

Genomic structural variants (SVs) are important sources of genetic diversity and contribute to the adaptation and evolution of plant-pathogenic fungi. Colletotrichum graminicola, the causal agent of maize (Zea mays L.) anthracnose, causes significant yield losses and serves as an important model for genetic studies. To investigate genomic plasticity in this pathogen, long-read genome assemblies were generated for four C. graminicola field strains collected in France, Croatia, Brazil, and Argentina and compared with the chromosome-level reference genome from the USA. Comparative analyses revealed extensive structural variation across core chromosomes, including large-scale translocations, inversions, duplications, and insertions and deletions, whereas minichromosomes showed markedly reduced collinearity among strains. Secreted proteins and predicted effectors were located significantly closer to regions affected by SVs. Approximately 174 kb of segmental duplications lacking genes were shared between core chromosome 04 and minichromosomes 11 and 13 in the reference genome but were absent from the other field strains, revealing a strain-specific structural link between core and minichromosomes. Protein clustering identified 733 candidate strain-specific genes across the reference genome. Approximately 72% of the genes annotated on the minichromosomes were strain-specific, indicating that most of their annotated genes lacked orthologous gene models in the other strains. Of the 733 candidate strain-specific genes, 492 retained homologous genomic sequences in at least one other strain and many also showed evidence of gene disruption. This pattern is consistent with lineage-specific pseudogenization and suggests that apparent gene presence-absence variation may arise through several evolutionary processes, including gene gain, gene loss and loss of gene function. Insertions and deletions in potentially regulatory flanking regions were associated with plant-pathogen interaction genes. In addition, interpretable machine learning analyses revealed a co-occurrence of duplications and LTR retrotransposons in genomic regions associated with secreted genes. This study shows that structural variation, transposable elements, and lineage-specific pseudogenization are important features of genome evolution in C. graminicola and identifies associations between genomic variation and genes involved in plant-pathogen interactions.

Approximately 8% of the human genome is derived from ancient retroviral infections. The most recently integrated of these endogenous retroviruses is the HERV-K(HML-2) clade, whose expression has been associated with cancer, amyotrophic lateral sclerosis, and embryogenesis. Studies of HERV expression, particularly HML-2, have relied predominantly on short-read sequencing. However, the high similarity among HML-2 proviruses prevents many short reads from being assigned uniquely to individual loci. We therefore compared haplotype-resolved long-read genome assemblies from 292 donors to resolve variation in proviral structure and coding capacity. Several loci previously thought to be fixed were structurally polymorphic. Tandem arrays occurred at 13 loci and contained up to six proviral copies in a single array. At 8q11.23, we identified a previously undescribed full-length provirus in one haplotype. All 583 other haplotypes carried a solo-LTR. We found that standard reference genomes failed to represent the coding capacity retained in many individuals, whose proviruses contained intact open reading frames despite disruptive mutations in the reference sequences. Short-read genotypes left 32.5% of the tested donor-variant pairs unresolved at sites associated with viral reading frames. These findings show why HML-2 expression must be interpreted in the context of the structural and coding alleles each individual carries.

Genetic rearrangement in the absence of selection is poorly understood. Long-read, high-depth sequencing has enabled analysis of the genetic landscape of spontaneous structural variants. We have applied PacBio Revio sequencing to identify the major classes of genetic alterations during long-term vegetative growth in Aspergillus flavus. These classes include a) precise deletions frequently associated with the flanking region of homology, b) clusters of short deletions that are heavily enriched within centromeric and subtelomeric regions that are typically 94% A+T, and c) deletions and expansions within a subset of homopolymer tracts. We find 70 homopolymer tracts longer than 44 bp, dispersed throughout the genome in noncoding regions and unassociated with transposable elements. A broader analysis of third-generation sequenced datasets from fungi found within the Phylum Ascomycota showed a high but variable prevalence of homopolymeric tracts across species. The adventitious identification of a putative flbA mutation, responsible for an altered conidial phenotype, by the Deep Variant AI system also demonstrates the power of this approach. These studies show that high-depth, long-read sequencing can probe SNPs, indels, homopolymeric tracts, and large-scale rearrangements in an unbiased fashion, and demonstrate that short deletions and homopolymer tract expansions arise spontaneously in A. flavus and other filamentous fungi.

Epstein-Barr virus (EBV) infects 90% of the population worldwide and is associated with multiple malignancies, including Burkitt, Hodgkin and post-transplant lymphomas, nasopharyngeal and gastric carcinoma. EBV also triggers autoimmune diseases, including multiple sclerosis and systemic lupus erythematosus. While host factors important for EBV B-cell entry have been defined, including CD21/CR2, MHC class II and R9AP, much remains to be learned about host factors that support post-entry steps of EBV infection. To gain insight, we conducted a human genome-wide CRISPR-Cas9 screen. In addition to known EBV B cell receptors/co-receptors, we identified multiple new factors and pathways essential for EBV B cell infection. These included multiple B cell signaling pathways, actin cytoskeleton and nuclear import factors. Knockout of CD19, its tetraspanin chaperone CD81 or karyopherin subunit alpha 1 (KPNA1, also called importin-5) significantly impaired establishment of EBV infection. EBV and CD19 co-localized at early stages of EBV infection, and the CD19 cytoplasmic tail was important for EBV infection, suggesting CD19 signaling supports EBV uptake. Depletion of KPNA1 significantly impaired establishment of EBV infection, but did not impair production of infectious EBV upon reactivation of latently infected B cells. Together, these results provide insights into pathways that support EBV B cell infection and identify potential therapeutic targets.

The Herpes Simplex Virus 1 (HSV-1) genome undergoes extensive chromatinization following entry into the nucleus, yet how viral regulatory elements influence this process during lytic infection remains poorly understood. The HSV-1 latency-associated transcript enhancer (LTE) is well established as a key element in the regulation of latency and reactivation, but its role during the productive infection has not been defined. Here, we identify a role for the LTE in regulating chromatin composition on the lytic HSV-1 genome. Deletion of a 310-bp region within the LTE selectively reduced enrichment of the canonical histone H3.1 at specific viral lytic genes, including ICP4 and ICP8, without altering H3.3 enrichment. This phenotype was not reproduced using a second LTE deletion that disrupted enhancer activity, indicating that regulation of H3.1 enrichment is separable from the function of the enhancer. LTE-dependent H3.1 enrichment persisted in the absence of viral DNA replication and was accompanied by reduced occupancy of the H3.1 histone chaperone CAF-1. However, depletion of CAF-1 did not substantially reduce H3.1 enrichment, indicating that additional mechanisms contribute to H3.1 association with the viral genome. Together, these findings reveal that the HSV-1 LTE regulates locus-specific chromatin composition during early lytic infection and identify a previously unrecognized connection between a viral enhancer and histone H3.1 enrichment on the incoming viral genome.

Archaeal extra-chromosomal elements (ECEs) play key roles in archaeal host biology and evolution, yet little is known about their diversity and interactions with each other. Here, we report 42 complete 66 - 165 kbp genomes: 23 "Klingons", 8 "Romulans" and "Romulan-like elements", and 11 "Ferengi", three new families of archaeal ECEs from diverse terrestrial environments. Methylation patterns, phylogenetics and gene sharing analyses suggest that Klingons and Romulans associate with Candidatus Methanoperedens and Ferengi are likely hosted by Candidatus Bathyarchaeaota. Like known archaeal viruses and plasmids, these ECEs both have circular genomes and replicate bi-directionally, unlike the linear Borg and mini-Borg ECEs that also replicate in Methanoperedens. The new ECE genomes all encode genes for methylation, glycosylation, protein degradation, and novel TnpB-like enzymes, which may interact with novel CRISPR-like arrays. All genomes also encode proteins with strong structural similarity to archaeal viral capsid proteins, and some are predicted to multimerize into pentameric or hexameric subunit capsid-like arrays and thus may encapsulate themselves. The three new ECE families share gene content with each other and with Borgs and mini-Borgs, often concentrated in specific genomic regions. Discordance among single gene and genomic phylogenies suggests that apparent relatedness of the five distinct ECE types stems from extensive horizontal gene transfer, rather than evolution from a common ancestor. The three new families of ECEs broaden our understanding of the diversity of archaeal ECEs and highlight the role of horizontal gene transfer in modulating the archaeal and archaeal ECE gene pools.

Lateral gene transfer (LGT) from organellar to eukaryotic genomes is ubiquitous, resulting in the presence of widespread nuclear mitochondrial DNA segments (NUMTs) and/or plastid gene transfers. While LGT from bacterial associates of eukaryotes is generally less common, many arthropod genomes are littered with LGT originating from Wolbachia (order Rickettsiales) - a vertically transmitted, obligate intracellular symbiont of invertebrates. This contrasts with Rickettsia, a related genus including many important human pathogens, which has not been shown to promulgate LGT within its arthropod vectors. Here, we present the 8.6 Gb genome of a cell line derived from A. variegatum (the tropical bont tick), which is the primary vector of Rickettsia africae (agent of African tick-bite fever). In addition to numerous NUMTs and diverse families of transposable elements, an almost-complete chromosome of R. africae origin was identified in the cell line genome, localised to the putative sex chromosome. Sequencing of field-collected A. variegatum confirmed the presence of this large LGT in wild ticks, although the absence of transferred genes from the R. africae plasmid provided a means to differentiate LGT from genuine rickettsial infections. These findings highlight the imperative to consider the possibility of LGT when screening vectors for pathogens by PCR.

Mobile genetic elements drive genome evolution and horizontal gene transfer, yet the mechanisms underlying the mobility of many eukaryotic elements remain poorly understood. Here we show that Starships, a recently discovered class of giant transposable elements in fungi, mobilise via an extrachromosomal circular intermediate. Using deep sequencing, we detect low-abundance reads spanning the hypothesised junction consistent with circularisation. We demonstrate that a synthetic circular intermediate, or "starmid", is sufficient for recombinase-dependent integration into the host genome and that circularisation enhances expression of the starship-encoded tyrosine recombinase through juxtaposition of 3' and 5' terminal sequences. This topology-dependent regulation parallels bacterial integrative elements, in which excision activates recombination machinery. Exploiting this circular intermediate, we explore starmids as a platform for transient DNA delivery, enabling CRISPR Cas9-mediated genome editing in Paecilomyces variotii with subsequent curing of the construct. Together, our results shed light on the biology of starship mobility and suggest promising biotechnology applications for these elements.

2A peptides are short (~18 to 22 amino acid) sequences utilized by a subset of RNA viruses to process their polyproteins in an unusual co-translational event known as StopGo or ribosome skipping. This process involves ribosome pausing at a conserved C-terminal D(V/I)ExNPGP motif, followed by peptidyl-tRNA hydrolysis of the nascent peptide-tRNAGly, and continued translation from the downstream proline codon. However, how nascent 2A peptide interacts with the ribosomal exit tunnel to promote these event(s) remains poorly understood and the molecular determinants governing 2A mediated StopGo translation remain undefined. Here, we developed a mammalian fluorescence resonance energy transfer (FRET) based high-throughput reporter system to quantitatively measure 2A activity across large libraries of natural and engineered 2A peptides. Systematic mutational analysis of canonical viral 2A peptides from porcine tescovirus and foot-and-mouth disease virus (P2A, F2A) identified key residues required for activity and revealed critical contributions from previously uncharacterized N-terminal residues. Extending this approach to metagenomic viral datasets uncovered conserved sequences governing StopGo efficiency leading to functional ranking of diverse viral 2A peptides in mammalian cells, including identification of a novel subclass that supports high stop-go activity (>99%). Phylogenetic analysis of viral genomes revealed clustering of functional 2A peptides and specific motifs within specific viral lineages branches, suggesting an evolutionary trajectory underlying 2A diversification and optimization. These findings define the sequence and functional landscape of viral 2A peptides and provide mechanistic and evolutionary insights into how nascent peptide ribosome exit tunnel interactions mediate StopGo translation.

Tandemly repeated sequences are noncoding, highly variable, yet integral components of eukaryotic genomes. Empirical research has so far largely focused on a few well-characterized, functional tandem repeats (TRs) or documented the broad diversity and rapid evolution of these repeats without explicitly testing the forces driving their evolution. By characterizing TR abundance and TR sequence composition within and between chromosomes and species in the stick insect genus Timema, we uncover that TR loads are enriched in genomic regions of high recombination but tend to be reduced in species with higher effective population sizes. These results suggest that there is a mutational effect of recombination that generates increased TR loads and that the bulk of TRs experience weak selection, leading to their rapid turnover over short evolutionary timescales. Finally, comparisons of centromere and non-centromere TRs reveal that the rapid divergence of centromeres is best explained by the combination of the fast-evolving nature of TRs per se and centromere-specific evolutionary forces. Overall, this study provides the first in-depth investigation of genome-wide TR evolution, while integrating the effects of recombination, genetic drift, and centromeric position along the genome.

Viruses generate immunogenic double-stranded RNA (dsRNA) during their life cycle, necessitating mechanisms to evade host innate immune detection. Protein kinase R (PKR) is a pattern recognition receptor (PRR) that senses viral dsRNA and activates the host integrated stress response (ISR) by phosphorylating eukaryotic initiation factor 2 (eIF2) on its alpha subunit, causing global mRNA translation arrest. Viruses evade activation of the ISR and translation shut-down to facilitate replication. Adenosine deaminase acting on RNA 1 (ADAR1) is uniquely positioned to support viral replication due to its ability to convert adenosine (A) into inosine (I) inside dsRNA, thereby disrupting base-pair matching and preventing PKR activation. Here, we discovered that Epstein-Barr virus (EBV), the oncogenic driver of nasopharyngeal carcinoma (NPC), requires ADAR1 both to maintain latent infection and for lytic reactivation by preventing PKR-mediated ISR activation. EBV lytic reactivation stimulated viral protein expression, which was attenuated in the absence of ADAR1. Importantly, lytic reactivation or treatment of synthetic dsRNA in the absence of ADAR1 induced phosphorylation of eIF2 by PKR, stress granule formation, and caused cancer cell lethality without the production of virion progeny. Together these findings support a strategy for EBV lytic induction therapy in NPC without viral dissemination by exposing EBV to endogenous PRRs through ADAR1 inhibition.